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SUMMARY 


Equations and procedures are presented for carrying out the numerical 
operations involved in computing the aerodynamic characteristics of lifting 
rotors by means of modern high-speed automatic computers. The character- 
istics considered herein include the thrust, profile drag, total power, 
flapping, rolling and pitching moments, direction of the resultant -force 
vector, and the harmonic contribution of the shear-force input to the 
hub. The equations are general and can account for stall and compressi- 
bility effects, hinge offset, and other factors that are normally omitted 
from conventional analytical rotor treatments. 


INTRODUCTION 


The problem of determining the aerodynamic characteristics of a 
lifting rotor basically consists of calculating the forces and moments 
on a rotating three-dimensional lifting surface from statically measured 
two-dimensional airfoil-section data. Because of the large number of 
degrees of freedom involved and the complications introduced by the vari- 
ation in section angles of attack and velocities over the rotor disk, an 
analytical approach to the problem necessarily contains various assumptions 
and simplifications. The largest number of such simplifying assumptions 
is found in the original efforts of Glauert and Lock, published almost 
30 years ago. Since that time, the theory has been progressively improved 
in accuracy and convenience of application by a number of investigators, 
and has been extended to meet the needs of new designs and more severe 
flight conditions . 

These extensions and improvements in accuracy, however, have been 
realized at the expense of a very large increase in the length and com- 
plexity of the resulting equations when closed-form solutions are desired. 
(See ref. 1, for example.) In order to make the application of these 
equations practical, it is usually necessary to evaluate the equations 
over a range of flight parameters and to present the results in chart 
form. (For example, see refs. 2 to 4 .) Such solutions, in spite of their 
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complexity, still do not account for blade stall outside of the reversed- 
velocity region, compressibility effects, and combinations of such design 
factors as hinge offset, blade twist and taper, and root cutout. 

Instead of trying to include such factors in the already complicated 
equations by first integrating the general equations and then substituting 
'articular values, it would seem more practical to account for these fac- 
ors numerically; that is, by substituting numerical values into the dif- 
■rential rotor equations and integrating numerically. This procedure is 
ie practicable only by the development in the last few years of high- 
■ed automatic computing machines, inasmuch as the numerical work involved 
evaluating a large number of cases by manual calculating methods is 
libitive. The availability of such computing machines to industry as 
as to research institutions, through ownership or rental, makes the 
•’ical procedures appropriate for design as well as for research 
Les . 

It is the purpose of this paper to present the equations and proce- 
ss necessary for numerically computing the aerodynamic characteristics 
Lifting rotors. The characteristics considered herein include the 
ust, profile-drag power, total power, flapping, rolling and pitching 
ents, direction of the resultant-force vector (necessary in stability 
.culations), and the harmonic contribution of each blade of the rotor 
the shear-force input to the hub . 


SYMB0K3 


q>B]_ coefficients of -cos \|r and -sin i|r, respectively, in expres- 
sion for 0; therefore, lateral and longitudinal cyclic pitch, 
respectively, deg 

i speed of sound, ft/sec 

a’ projection of angle between rotor resultant -force vector and 

shaft axis in plane containing flight path and shaft axis, 

tan“l — 

Cip 

&q constant term in Fourier series that egresses p, radians; 

hence, rotor coning angle 

a^ coefficient of cos n\|r in expression for p 

B tip-loss factor; blade elements outboard of radius BR are 

assumed to have profile drag but no lift 
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b' 


*n 

C 0 

C H 


- , m 


'n 


°Q 

Cip 

Crj(i|r ) 


°Y 


number of blades per rotor 

projection of angle between rotor resultant -force vector and 
shaft axis in plane perpendicular to plane containing flight 

1 Cy 

path and shaft axis, tan _J - — 

c T 

coefficient of sin n\|r in expression for (3 
constant term in Fourier series that expresses Mr I h a2 

H 


rotor longitudinal-force coefficient, 


jtR 2 p(0R) : 


rolling -moment coefficient. 


pitching-moment coefficient. 


L' 


nR 2 p(fiR) 2 R 


M 


jtR^nR)^ 

coefficient of cos n\|r in expression for 

P 


rotor-shaft power coefficient. 


rotor-shaft torque coefficient. 


nR p(ar) : 
Q 


rtR 2 p(ftR) 2 R 


rotor thrust coefficient, 


T 


nR 2 p(ftR) 2 


Crp value at a given azimuth position \|r, expressed by series 
C T ( t ) = E 0 + E X cos \|r + F-|_ sin \jr + E 2 cos 2\(r + F 2 sin 2i|r + 

cos + F^ sin 

shear-force coefficient at flapping hinge due to mass forces, 

F m 


2 2 
nR p(ffi) 


rotor lateral-force coefficient. 


nR 2 p(aR ) 2 

blade-section chord at radial station x, ft 
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section profile-drag coefficient 


equivalent blade chord (on thrust basis). 


■'BR 


cr^ dr 


■'BR 


ft 


r^ dr 


equivalent blade chord (on thrust-moment basis). 


n BR 

/ cr^ dr 


-'BR 


-, ft 


c-5 dr 


''l 

c t 

»n 

E 0 

E n 


m 


•n 


H 


L h 


k = 


section lift coefficient 
blade chord at tip, ft 

coefficient of sin n\|r in expression for 

constant term in series expression for Crp(i|r); = Cp 

coefficient of cos n\|r in expression for Cp(\|r) 

offset of center line of flapping hinge from center line of 
rotor shaft, ft 

shear force at flapping hinge (perpendicular to blade-span 
axis) due to mass forces, lb 

coefficient of sin ni|r in expression for Cp(i|r) 
gravitational acceleration, ft/sec^ 

component of rotor resultant force perpendicular to rotor shaft 
in longitudinal plane of symmetry, lb 

mass moment of inertia of blade about flapping hinge, 

m(r - e)^ dr, slug-ft^ 


em b 
Mm, s 


L’ 


rotor rolling moment, lb -ft 
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M 

^n,s 

Miji 

% 

M x 

m 

mb 

P 

P 

Q 

q. 

R 

r 

T 

u 

Up 

"Ujji 

v 

v 

x 

Y 

a 


rotor pitching moment, lb -ft 

mass moment of blade about center line of rotor shaft, slug-ft 
thrust moment of blade about flapping hinge, lb-ft 
weight moment of blade about flapping hinge, 


Mach number of blade section, ufffi/a 
mass of blade per foot of radius, slugs/ft 
mass of blade, slugs 
shaft power, ft-lb/sec 

helicopter rolling velocity, radians/sec 
shaft torque, lb-ft 

helicopter pitching velocity, radians/sec 

blade tip radius measured from center of rotation, ft 

distance measured along blade from axis of rotation to blade 
element, ft 

component of rotor resultant force acting along rotor shaft, lb 

nondimens ional resultant velocity perpendicular to blade-span 
axis at blade element 

component at blade element of u parallel to shaft axis 

component at blade element of u perpendicular to blade-span 
axis and to shaft axis 

true airspeed of helicopter along flight path, ft /sec 

induced inflow velocity at rotor, ft/sec 

ratio of blade-element radius to rotor-blade radius, r/R 

lateral component of rotor resultant force perpendicular to 
both T and H, lb 

angle between shaft axis and plane perpendicular to flight 
path, positive when axis is pointing rearward, deg 
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Or 

■blade-element angle of attack, measured from line of zero lift, 
deg 

P 

blade flapping angle with respect to shaft at a particular 

azijnuth position, - a p cos V ~ ^l s ^ n ^ “ a 2 cos ^\|r • • • > 

radians 

7 ' 

PCp 'R^ 

mass constant of blade, — 

h. 

7 ' 

' X 

1|_ 

local value of blade mass constant, 

J-h •> 

P M w 
TJ = £. 

s X h 


0o 

collective-pitch angle at blade root, average value of instanta- 
neous blade-root pitch angle around azimuth, deg 

01 

difference between blade-root and blade-tip pitch angles, posi- 
tive when tip angle is larger, deg 

9 .75 

blade-section pitch angle at O .75 radius, deg 

0 

instantaneous blade-section pitch angle; angle between line of 
zero lift of blade section and plane perpendicular to rotor 
shaft, 0 q + 0px - A]_ cos i|r - Bj_ sin \|r, deg 

A 

inflow ratio, - si ^ R a ~ ~ 

M- 

tip-speed ratio, - ^° s a 

1 = e/R 


P 

mass density of air, slugs/cu ft 

a 

rotor solidity, bc e /jtR 

a x 

1 

local solidity, bc/rtR 

4 

inflow angle at blade element in plane perpendicular to blade- 

1 3 - A 

span axis, tan ^ — , deg 

Uip 


span 


A 
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blade azimuth angle measured from downwind position in direction 
of rotation,, deg 

& rotor angular velocity, radians/sec 

Subscripts: 

c radius of cutout; that is, radius at which lifting surface of 

blade begins 

a accelerating 

d decelerating 

0 profile 

1 induced 

l lift 


METHOD OF ANALYSIS 


The process of numerically determining the aerodynamic character- 
istics of a rotor operating at a particular flight condition consists of 
calculating the individual force (and moment, if desired) contributions 
of a specific number of blade sections at various points on the rotor 
disk, averaging the values around the disk at a particular radial station, 
and then radially integrating these averages along the blade to obtain 
the result. The equations required to evaluate the force contribution 
of a blade section and the procedure for utilizing the equations are 
written specifically for a hinged (i.e., free-to-cone ) rotor, but can be 
applied with but little modification to seesaw rotors. The modification 
of the equations to the seesaw-rotor configuration is discussed in refer- 
ence 5- 


Reference-Axis System 

In the analysis developed herein, the flapping hinge may be offset 
radially from the center line of the shaft. For such a configuration, 
flapping and feathering motions are not directly equivalent and it is 
therefore most convenient to reference all angles and velocities to an 
axis coinciding with the center line of the rotor shaft. The rotor shaft 
is therefore used as the reference axis in this paper instead of the no- 
feathering axis which was used in references 1 to 4. 
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Equations 


t 


Most of the equations listed in this section are based on, or can he 
readily derived from, standard blade-element equations. (See refs. 1 
and 6, for example.) Because the equations are set up for numerical inte- 
gration, most of the aerodynamic restrictions that are normally employed 
in rotor analyses are avoided. Thus, the equations can be used for blades 
of any airfoil section, mass distribution, twist, plan-form taper, root 
cutout, and flapping-hinge offset and can account for the effects of heli- 
copter pitching and rolling velocities. Flapping -hinge configurations 
wherein blade flapping and lagging motions result in pitch-angle changes 
can be simply considered by the addition of the appropriate terms in the 
equation for section angle of attack presented in this section. Stall 
and compressibility effects can be accounted for to the extent of using 
actual airfoil data that correspond to the range of angle of attack and 
Mach number encountered at all points of the disk, including the reversed- 
velocity region. Although no small-angle limitations have been made 
regarding the section inflow angles or the treatment of the magnitude of 
the resultant velocities and the direction of the resultant— force vectors 
at the blade sections, the analysis does incorporate small-angle assump- 
tions in regards to the blade flapping angles. In addition, it is assumed 
that the rotor blades are rigid in bending and torsion and that radial- 
flow effects can be neglected. The equations can also be applied to any 
known or assumed radial or azimuth variation of induced velocity. 

Input quantities .- The input parameters fixing the flight condition 
are A, 0 anri p. These parameters may be either assumed directly 

as independent variables or estimated from known parameters (such as 
thrust and shaft power) by means of published rotor-performance charts 
(such as those of refs. 2 and 3). In addition to the flight parameters, 
the following design and operating variables are assumed to be known: 

0 1 , A-l, Bj_, c(x), I h , e, M w , n, b, R, p, p, q, and m(x), as 

well as two-dimensional blade-section airfoil data. 

Equations for the calculation of section velocities, angles of attack, 
and Mach number.- Equations that are necessary for the calculation of sec- 
tion velocities, angles of attack, and Mach number are listed as follows: 


u T = x + p sin \|r 


( 1 ) 



(x - |)(2a 2 sin 2\|r - 2b 2 cos 2i | r + 3a^ sin 3t - 3bj cos 3t) 
p cos \|r(a 0 - a 1 cos xir — b^ sin \|r - a g cos 2\|r - b 2 sin 2\|r - 
a^ cos 3t “ *>3 sin 3ty) 


( 2 ) 
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OCj, = 0 + 0 


u 2 = Urp 2 + Up 2 

(3) 

0 = tan” 1 ^ 

U-rp 

(4) 

x - Ap cos \|r - Bp sin \|r + 0 

(5) 


m y = 


_ uflR 


a 


( 6 ) 


Equations for the calculation of flapping coefficients .- The flapping 
coefficients may be calculated from the following relations: 


M, 


Ih°' 


T _ 1 

5 2 




7' u 2 (x - |) 


c £ cos 0 & x + 


. 1.0 


x r 


/' x u 2 (x - l)c d ^ 0 sin 0 dx 


(7) 


% 


i h a 


■g = Cq + cos \(r -H D-^ sin \|r + Cg cos 2\|r + s ^ n 2\|r + cos 3^ + 

D 5 sin 3* (8) 


.L.JiL 


a - * ' c o - r 1 *- ' — 

\ J -h 


'°' l ° I ^ 2 I 1+ " 


(9) 


Aa-^ = 


_ D i + [jlfbi + £bi) - 2 |(l + r)g)J 


(n£l) (lOa) 


Aa n = 


P 1 + [^( b l + ^l) - 2 + Tlg)j 


(h > 1) (10b) 


T 7 a 
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z^b 1 = 


C 1 + [jl( a l + ^l) + 2 §(! + ng)] 


SI 

8 


(b 2 + | v^jy' 


a 


(id 


a~ = 


3 - ti 
D 2 

3 - T] 
C 3 

8 - T) 

D 3 

3 8 - r, 


t>« = 


a 2 = 


= 


( 12 ) 

(13) 

(14) 

(15) 


Equations for the calculation of performance parameters . - The per- 
formance parameters may he determined from the following equations: 


s?), * 5 “Ai 008 ** 


(16) 


JF L ■ 5 'A,, •*» * 


(17) 


' 1,1 


' B L ^ /dC 


n 


x c \|r=l 


'T 


dx )i 


dx 


(18) 


Pl.o *50. / dC \ 


(19) 


fcS dCm 

where refers to the summation of — +■ values computed at n azimuth 

♦=1 


dx 


4. X- (-B I , A 0 360° 2 ( 360 O ) 

stations. (For example, \|r = 0 , ^ , 


n 


n 


•) 
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Also 


C T = C T,2 + C T,o 

(20) 

— = i a u^xc 7 sin 0 
dx 2 x l * 

(21) 

— = — a x u^xc d n cos 0 
dx 2 x d ^° 

(22) 

PB t=n an 

n - / IV — ^ dx 

u Q,a J n 2_ dx 

x c t=l 

(23) 

pl’O , dC Q d 

c Q,d - / „ X -df 1 - 

Ox c i|r=l 

(24) 

C Q = C Q,d “ C Q,a 

(25) 

- 1 a u 5 c, 
dx 2 x d,o 

(26) 


1.0 , t=n 


j p j. . V T 

X \|r=l 


dC P,o 


dx 


dx 


(27) 


Equations for the calculation of rotor forces and moments .- The fol- 
lowing equations may be used to calculate certain rotor force and moment 
coefficients that may be useful in vibration and stability analyses: 


dCg ± ~ 2 °x u ^ c i(~ s ^ n $ s: *- n ^ “ cos 0 sin 3 cos 


(28) 


dCg^ 0 = i a x u^c^ Q (cos 0 sin \|r - sin 0 sin 3 cos \|f)dx 


(29) 


, . r B i dc H,i 


(30) 
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,1.0 \|f=n fln 

1 V- d U H,o 


n 


'x r 


t=l 


dx 


dx 


C H “ C H,i + c H,o 


dCy i = ^ a x u 2 c^(sin 0 cos \|r - cos 0 sin p sin \lr)dx (33) 


dCy o = \ cr x u^C£ o^“ cos ^ cos ^ “ s -*- n ^ s ^ n P s ^ n A l r ) < ^ x (3^) 


~L “ 


B \|r=n 


dC Y,i 


t=l 


dx 


dx 


,1.0 \|r=n 


°Y,o - 


dC Y,o 


/ n dx 

x c f 1 


dx 


C Y “ C Y,i + C Y,c 


a' = tan“-*- 


b ' = tan -1 - ^ 


(31) 

(32) 


(35) 

(36) 

(37) 

(38) 

(39) 


C T (i|r) = a x u2c Z cos 0 dx + J a x u2c d^o sin 0 (^°) 

\ x c x c ) 

After Ciji(,|r) is evaluated at different azimuth stations, a harmonic 
analysis would express Cgi(\[r) as a function of \|r as follows: 

C T (i|r) = E Q + cos \|r + F-^ sin \|r + Eg cos 2i|f + Fg sin 2\(r + 

E^ cos 3\|r + F^ sin 3t (4l) 
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Then 




\ 


£z = _ 

a 2a 


(42) 


Qm _ _ ^1 
0 ~ 2o 


(43) 


Equations (42) and (43) represent aerodynamic rolling and pitching 
moments. Blades with offset flapping hinges also produce hub moments 
due to centrifugal forces as given by the following equations: 



JL 1 + j 
7' 1 - | 


(44) 



1 + I 
l - I 


(45) 


The oscillating shear force at the flapping hinge (perpendicular to 
the blade-span axis) is composed of thrust and mass forces. The thnjst- 
force contribution, expressed in harmonic form, is given by equation (4l). 
The nondimens ional mass-force contribution, made up of centrifugal and 
inertia forces, is represented by 


"bMjji c* r~~ 

Cy = - — j^- (l) 2 — j^a^ - a^ cos \|r - b^ sin \|r - ag cos 2ijr - \>2 sin 2\|r - 
a 3 cos 3* - b^ sin 3t) + (l - k)(a x cos y + b x sin y + 

4a2 cos 2 \| r + 4bg sin 2\|r + 9a^ cos 3\Jr + 9b^ sin 3\|r]j 


(46) 


Inasmuch as equations (4l) and (46) contain no aerodynamic damping due to 
bending or bending inertia terms, they apply only to inflexible blades. 


Calculation Procedure 

The procedure for calculating rotor characteristics is outlined as 
follows: 

(l) For the given input quantities, assume or calculate an approxi- 
mate set of flapping coefficients. (The charts of ref. 4 may be used to 
calculate these approximate coefficients if a^ and b^ are assumed to 

be equal to zero. ) 


% 


Ik 
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( 2 ) With the given input quantities and the approximate flapping 
coefficients , calculate section velocities and angles of attack by means 
of equations (l) to (5). 

( 5 ) Calculate section Mach numbers by means of equation (6). 

(t ) By using the calculated angles of attack and Mach numbers, 
obtain the corresponding values of c^ and c^ Q from appropriate two- 

dimensional section data. 

(5) Calculate Mij /l h fi 2 by means of equation (7), make a harmonic 
analysis of the result, and express the answer in the form of equation (8). 

(6) Calculate new values of flapping coefficients by means of equa- 
tions (9) to (15). (The incremental first-harmonic coefficients given by 
eqs. (10) and (ll) should be added as correction factors to the output 
values from the preceding iteration.) 

(7) Repeat steps (2) to (6) until the output values of flapping 
coefficients differ negligibly from the input values. 

(8) Compute 0 for each value of \| r, using the final values of 
flapping coefficients. 

(9) With the use of the lift and drag coefficients and velocities 
corresponding to the final set of flapping coefficients, calculate the 
performance and stability parameters by means of equations (l6) to ( 45 )- 

It will be noted that additional performance, stability, and loads 
parameters can be calculated with little additional effort, inasmuch as 
the section velocities, angles of attack, and lift and drag coefficients 
have already been computed. 


DISCUSSION OF METHOD 
Assumptions of Theory 


Although the numerical integration procedure eliminates many of the 
restrictive assumptions and simplifications of conventional rotor theory, 
a number of assumptions still remain. The degree to which these suppo- 
sitions can affect the validity of the results depends on the specific 
applications of the method. Perhaps the most important assumption, and 
one which is inherent in all the theoretical rotor treatments developed 
thus far, is the use of steady, two-dimensional airfoil-section character- 
istics in evaluating rotor forces and moments . Although such an approach 
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has been proven successful in predicting rotor characteristics under 
normal operating conditions, a complete evaluation of its accuracy in 
\ conditions involving large amounts of stall and compressibility effects 

must await further experimental checks . Limited theory-data comparisons 
that have been made thus far for extreme operating conditions, however, 
do show good agreement. The application of such theory to these condi- 
tions is therefore felt to be the best step which can be taken for pre- 
liminary estimates at the present time. As more experimental data become 
available, comparison with these predictions will serve to show whether 
or not empirical or other revised approaches are necessary. 

The assumption of blade inflexibility in bending and torsion would 
„ normally be unimportant except in aeroelastic problems. Bending flexi- 

bility should be considered, for example, in the computation of coupled 
blade -fuselage vibrations, whereas torsion flexibility would enter into 
\ blade -fuselage clearance calculations or in stability or stress calcula- 

tions wherein a large blade-section center-of -pres sure shift associated 
with operation at high Mach numbers occurs. The neglect of radial-velocity 
components along the blade would probably be significant only at very high 
tip-speed ratios wherein they might influence the profile-drag power and 
the stall pattern of the rotor. 

An important limitation inherent in the analysis, and one whose 
effects can be readily evaluated, is the use of small-angle assumptions 
as applied to the rotor-blade flapping angle. The flapping analysis 
containing these assumptions yields erroneous answers if the amplitude 
of the flapping angle is in the neighborhood of 15° or above. For flight 
conditions involving larger amplitudes, the numerical method for calcu- 
lating blade flapping presented in reference 5, which contains no small- 
angle limitations, should be used. (The two methods yield identical 
answers for conditions wherein the flapping angles are below approxi- 
mately 15° • ) The method of reference 5 can also be used for cases in 
which a knowledge of the transient blade-flapping behavior in maneuvers 
is desired. 


Application of Method 

1 

Aside from the assumptions just discussed, the accuracy of the 
answers obtained from the method of this paper is dependent upon the 
number of radial and azimuth stations around the disk used in the inte- 
gral force, moment, and power equations. The problem is to determine 
the point of diminishing returns at which further increases in accuracy 
obtained through the use of a greater number of stations are counter- 
balanced by the additional time and cost required to compute the extra 
stations. The answer to this problem is a variable and is dependent upon 
the flight condition for which the answer is to be calculated (i.e., the 
extent of the regions on the disk undergoing stall, compressibility, and 
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reversed -velocity effects), the type of answer to be calculated (for 
example, thrust loading for bending -moment calculations, stability deriv- 
atives, or blade flapping), and the type of automatic computing equipment 
available for the calculations. 

For example, because the calculation of the flapping coefficients, 
which is necessarily an iterative procedure, is quite lengthy, it is 
important to determine the initial estimate as closely as possible and 
to require no more accuracy than is necessary for the particular appli- 
cation. Thus, for cases in which the determination of the flapping coef- 
ficients is not in itself the objective but is merely the means for calcu- 
lating such items as rotor thrust or power, it is usually sufficiently 
accurate to ignore the third, and perhaps even the second, harmonics of 
flapping, and to keep the number of iterations at a minimum. Such con- 
siderations become less important, of course, when the calculations are 
carried out on a modern high-speed automatic computer. Conversely, if the 
flapping motion is desired in its own right, the number of iterations 
(and the number of harmonics) can be carried out to any desired accuracy. 

In order to obtain an idea of the effect of the number of stations 
on the accuracy of the computations, rotor flapping coefficients (called 
part I solutions) and thrust, power, and direction of the resultant-force 
vectors (called part II solutions) were evaluated for a sample flight 
condition by using different numbers of radial and azimuth stations. 
Angle-of-attack and Mach number contours corresponding to the sample case 
are plotted in figure 1 to illustrate the rate of change of these quan- 
tities with radius and azimuth angle. As can be seen from the figure, 
the flight condition chosen was rather extreme and involved large amounts 
of stall and compressibility effects at a tip-speed ratio of 0.50 and a 
tip speed equal to 75O feet per second. The results of the calculations 
are presented in table I. 

The results shown in table I were obtained at thirty-six (Aa(r = 10°), 
eighteen (Ai|r = 20°), and eight (A\|r = 45°) azimuth stations. In the radial 
direction, values were computed at x = 0.15 and at x = 1.0, with seven 
(Ax = O.IO 63 ) and three (Ax = 0.2125) uniformly spaced stations in between. 
(For those quantities, such as thrust and induced torque, which require 
the integrations to be carried out to x = 0.97 instead of x = 1.0, 
values were also computed at x = 0.97- The value of the quantity between 
x = 0.97 and x = 1.0, computed by the trapezoidal rule, was subtracted 
from the total value of the quantity obtained by integrating from x = 0 . 15 
to x = 1.0.) With the exception noted in the table, all the part II 
solutions were computed with the same set of flapping coefficients (those 
corresponding to Ax = 0.2125, A\|r = 45°) in order to isolate the effect 
of changes in the number of stations on the results . 

It can be seen that, even for the extreme flight condition considered, 
the differences in flapping angles resulting from the use of a different 
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number of stations are small, being of the order of a few tenths of a 
degree for the coning angle and first harmonics, and even less for the 
second and third harmonics. (For these cases, the number of iterations 
was such that the values listed in the table differed from those of the 
preceding iteration by less than 0.0 6° for aQ, a^, and b-j_ and less 

than 0.03° for b2> and 

Table I shows that the part II solutions are somewhat more sensitive 
to the number of stations used than are the flapping coefficients, although 
not unduly so. By using the results for the case wherein Ax = O.IO63 
and = 10° as a base, it can be seen that approximately doubling the 
interval of radial and azimuth stations results in an error of only 1 per- 
cent or less in the thrust and power values. When A\|r is increased from 
20° to 4-5°, the error jumps to the order of 5 percent , which, in many cases 
(such as those involving extensive stall and compressibility effects), is 
within the order of accuracy and use of the airfoil-section data. In cer- 
tain calculations, however, such as comparative studies, the introduction 
of a 5 “Per cent error might be undesirable and could be avoided by doubling 
the number of azimuth stations. In general, it appears that greater 
accuracy is achieved with a given increase in stations by increasing the 
number of azimuth, rather than radial, stations. 

In order to illustrate the effect of differences in flapping coeffi- 
cients on the part II solutions, the case for Ax = 0.2125 and Aty » 4-5° 
was repeated, but with the second and third flapping coefficients assumed 
equal to zero. As can be seen in table I, there is a slight increase 
(about 2 percent) in thrust, and less than one-half of 1 percent decrease 
in power. 

In summary, it appears that if accuracy is a prime consideration, 
or if a very-high-speed computer (such as an IBM type 704- electronic data 
processing machine or a Remington Rand UNIVAC) is available, then approx- 
imately five radial and eighteen azimuth stations should be used for most 
part II solutions, whereas the part I (flapping) solutions can be deter- 
mined with the same accuracy by using five radial and only eight azimuth 
stations. For slower speed computers, such as the IBM Card-Programmed 
Electronic Calculator, five radial and eight azimuth stations can be used 
for the part II solutions with a considerable saving in computing time 
at the expense of a small reduction in accuracy. 

For a given number of stations, the possibility exists that the 
accuracy of a set of calculations can be improved by nonuniform spacing 
of the radial and azimuth stations. The possibility is greatest in 
regards to the radial spacing, particularly when a small number of sta- 
tions, such as five, is utilized. The effectiveness of a specific choice 
of the number and spacing of radial stations is also dependent on the 
manner in which the radial integrations are carried out. Simpson's rule 
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was used in the sample calculations presented in table I and was found 

to be satisfactory for the type of thrust and power loading curves 

encountered. / 

The technical procedures involved in setting up calculations of the 
type presented in this paper are described in a number of references 
(see ref. 7 > for example) and will not be discussed herein. It is con- 
sidered worthwhile to point out, however, that approximately l60 hours 
and 120 hours of programming (i.e., problem setup) time are required for 
the part I and II solutions, respectively. This part of the work is done 
only once, of course, and need not be repeated for each new set of input 
variables. Approximately 18,000 to 20,000 operations are involved in the 
part I solutions for the 48-station case when compressible section data 
are included, and approximately an equal number for the part II solutions. 

The corresponding time required to run off a part I iteration or a part II 
case for each new set of variables would vary from approximately 100 minutes f 
on an IBM CPC to approximately 20 seconds on an IBM type 704 calculator. 

The machine time required per case varies almost directly with the number 
of stations used. 


CONCLUDING REMARKS 


A numerical method suitable for application to automatic computing 
machines has been presented for calculating the aerodynamic character- 
istics of rotors for flight and design conditions that are outside the 
range of conventional analytical rotor theory. By means of the equations 
presented, the effects of such items as stall, compressibility, flapping- 
hinge offset, blade mass factor, twist, taper, and root cutout on such 
rotor characteristics as flapping, thrust, power, and hub moments, as well 
as on certain stability derivatives, can be investigated. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 6, 1956. 
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TABLE I.- EFFECT OF NUMBER OF RADIAL AND AZIMUTH STATIONS ON CALCULATED ROTOR 
CHARACTERISTICS FOR A SAMPLE FLIGHT CONDITION 
! vj. = O.5O; A = -0.091; 9 0 = 15°; 0 L = -8°; GR = 75O ft/sec; 7' = 1.0; £ = 0; x c = O.I5] 


(a) Part I solutions 


Ax 


^Q> 

a l> 

b l> 

& 2 , 

b 2 , 

a. 3 . 

b 3> 

deg 

deg 

deg 

deg 

deg 

deg 

deg 

deg 


10 

2.23 

11.34 

2.67 

1.022 

0.114 

-0 . 001 

0.186 

0 . 1063 

20 

2.22 

11.32 

2.70 

1.017 

.117 

-.001 

.184 


45 

2.14 

11.16 

2.83 

• 938 

.179 

.056 

.144 


10 

2.18 

10.91 

2.63 

1.021 

.121 

.004 

.188 

.2125 

20 

2.18 

IO.96 

2.65 

1.017 

.125 

-.001 

.186 

45 

2.18 

11.21 

2.79 

•959 

.183 

.058 

• 153 


(t>) Part II solutions 8. 


Ax 

A\(r, 

deg 

Cip/CT 

ACt y 
percent 

c P,o l a ’ 

AC P,o 1°’ 
percent 

Cp/a 

ACp jo , 
percent 

a ' , 

deg 

Aa' , 
deg 

deg 

Ab ' , 
deg 


10 

0.0685 

______ 

O.2151 

_____ 

0.1761 

— 

13-75 

— 

-0.47 

— 

0.1063 

20 

.0697 

1.8 

.2130 

1.0 

.1785 

1.4 

13.99 

0.24 

-.32 

0.15 

45 

.0654 

- 4.5 

.2278 

5.9 

.1842 

4.6 

14.39 

.64 

-1.08 

-.61 


10 

.0673 

-1-7 

.2169 

.8 

.1745 

-.9 

13.69 

-.06 

-.52 

-.07 

.2125 

20 

.0686 

.1 

.2151 

0 

•1775 

.8 

13.94 

.19 

-.36 

.11 

45 

.0660 

-3-7 

.2289 

6.4 

.i860 

5.6 

14.40 

.65 

-1.16 

-.69 


45 

.0672 

-1-7 

.2158 

• 3 

.1765 

.2 

14.11 

.36 

.02 

.49 


a All part II solutions (except the last one) were calculated with flapping coefficients 
corresponding to Ax — 0.2125 and A\|/ = 45 1 ” 1 . The last solution was calculated with same 
values of aQ, ap, and b]_, but with a 2 , b 2 , a 3> b 3 assumed equal to zero. 
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Drag divergence Mach 
number exceeded 

Stall angle exceeded 


Figure 1.- Angle of attack and Mach number contours for sample flight 
condition, n = 0.50; A = -0.091; e Q = 35°$ e i = " 8 °» = 75° feet 

per second; 7 ' = 1 . 0 ; £ = 0 ; x c = 0 . 15 . 


NACA - Langley Field, Va. 


